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Abstract 

We present a systematic study on linear propagation of ultrashort laser 
pulses in media with dispersion, dispersionless media and vacuum. The 
apphed method of amplitude envelopes gives the opportunity to estimate 
the limits of slowly warring amplitude approximation and to describe an 
amplitude integro-differential equation, governing the propagation of opti- 
cal pulses in single cycle regime. The well known slowly varying amplitude 
equation and the amplitude equation for vacuum are written in dimension- 
less form. Three parameters are obtained defining different linear regimes 
of the optical pulses evolution. In contrast to previous studies we demon- 
strate that in femtosecond region the nonparaxial terms are not small 
and can dominate over transverse Laplacian. The normalized amplitude 
nonparaxial equations are solved using the method of Fourier transforms. 
Fundamental solutions with spectral kernels different from Fresnel one are 
found. One unexpected new result is the relative stability of light pulses 
with spherical and spheroidal spatial form, when we compare their trans- 
verse enlargement with the paraxial diffraction of lights beam in air. It 
is important to emphasize here the case of light disks, i.e. pulses whose 
longitudinal size is small with respect to the transverse one, which in some 
partial cases are practically diffractionless over distances of thousand kilo- 
meters. A new formula which calculates the diffraction length of optical 
pulses is suggested. 

1 Introduction 

For long time few picosecond or femtosecond (fs) optical pulses with approxi- 
mately equal duration in the x, y and z directions (Light Bullets or LB), and fs 
optical pulses with relatively large transverse and small longitudinal size (Light 
Disks or LD) are used in the experiments. The evolution of so generated LB and 
LD in linear or nonlinear regime is quite different from the propagation of light 



beams and they have drawn the researchers' attention with their unexpected 
dynamical behavior. For example, self-channeling of femtosecond pulses with 
power little above the critical for self-focusing [T] and also below the nonlinear 
collapse threshold 2 (linear regime) in air, was observed. This is in contra- 
diction with the well known self-focusing and diffraction of an optical beam 
in the frame of paraxial optics. Various unidirectional propagation equations 
have been suggested to be found stable pulse propagation mainly in nonlinear 
regime (see e.g. Moloney and Kolesik [3], Couairon and Mysyrowicz [4], Chin at 
all. [5], for a review). The basic studies in this field started with the so called 
spatio-temporal nonlinear Schrodinger equation (NSE) which is one compila- 
tion between paraxial approximation, the group velocity dispersion (GVD) and 
nonlinearity O [71 HI |9] . The influence of additional physical effects were studied 
by adding different terms to this scalar model as small nonparaxiality [HI [18] , 
plasma defocussing, multiphoton ionization and vectorial generalizations. It is 
not hard to see that for pulses with low intensity (linear regime) in air and gases 
the additional terms as GVD and others become small and the basic model can 
be reduced to paraxial equation. This is the reason diffraction of a low intensity 
optical pulse governed by this model on several diffraction length to be equal to 
diffraction of a laser beam. On other hand, the experimentalists have discussed 
for a long time that in their measurements the diffraction length of an optical 
pulse is not equal to this of a laser beam z^^j™ = k^r^, even when additional 
phase effects of lens and other optical devices can be reduced. Here ko denotes 
laser wave-number and r'j_ denotes the beam waist. Thus exist one deep differ- 
ence between the existing models in linear regime, predicting paraxial behavior 
in gases, and the real experiments. 

The purpose of this work is to perform a systematic study of linear propaga- 
tion of ultrashort optical pulses in media with dispersion, dispersionless media 
and vacuum and to suggest a model which is more close to the experimental 
results. In addition, there are several particular problems under consideration 
in this paper. 

The first one is to obtain (not slowly varying) amplitude envelope equation 
in media with dispersion governing the evolution of optical pulses in single-cycle 
regime. This problem is natural in femtosecond region where the optical period 
of a pulse is of order 2 — 3 fsec. The earliest model for pulses in single-cycle 
regime suggested by Brabec and Krausz l20j is obtained after Taylor expansion 
of wave vector k'^{uj) about ujq. It is easy to show [21] that this expansion 
diverge in solids for single-cycle pulses. The higher order dispersion terms start 
to dominate and the series can not cut off. This is the reason more carefully 
and accurately to derive the envelope equation before using Taylor series. In 
this way we obtain an integro-differential envelope equation where no Taylor 
expansion of the wave vector k'^{uj), governing evolution of single cycle pulses 
in solids. 

The second problem is to investigate more precisely the slowly varying en- 
velope equations governing the evolution of optical pulses with high number of 
harmonics under the envelope. The slowly varying scalar Nonlinear Envelope 
Equation (NEE) is derived in many books and papers [TT ] [T2 ] [TS ] [M ] [T5l [T6[ [T7] . 
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After the deriving of the NEE, most of the authors use a standard procedure 
to neglect the nonparaxial terms as small ones. Only some partial nonparaxial 
approximations in free space [inillllllH] and optical fibers [H] were studied. In 
P2] we rewrite the NEE in dimensionless form and estimate the influence of the 
different linear and nonlinear terms on the evolution of optical pulses. We found 
that both nonparaxial terms in NEE, second derivative in propagation direction 
and second derivative in time with coefficient, are not small corrections. 
In fs region they are of same order as transverse Laplacian or start to dominate. 
These equations with (not small) nonparaxial terms are solved in linear regime 
|22| and investigated numerically in nonlinear [23j . In this paper we include 
GVD term in the nonparaxial model and study also the envelope equation of 
electrical field in vacuum and dispersionless media. It is important to note that 
the Vacuum Linear Amplitude Equation (VLAE) is obtained without any ex- 
pansion of the wave vector. That is why it work also for pulses in single-cycle 
regime (subfemto and attosecond pulses). 

Last but not least the nonparaxial equations for media with dispersion, dis- 
persionless media and vacuum are solved in linear regime and new fundamental 
solutions, including the GVD, are found. The solutions of these equations pre- 
dict new diffraction length for optical pulses z^f^f = kgr'^/zQ, where Zq is the 
longitudinal spatial size of the pulse (the spatial analog of the time duration 
to; Zq = vto; V is group velocity). In case of fs propagation in gases and vac- 
uum we demonstrate by these analytical and numerical solutions a significant 
decreasing of the diffraction enlargement in respect to paraxial beam model and 
a possibility to reach practically diffraction-free regime. 

2 From Maxwell's equations of a source-free, dis- 
persive, 

nonlinear Kerr type medium to the amplitude 
equation 

The propagation of ultra-short laser pulses in isotropic media, can be charac- 
terized by the following dependence of the polarization of first Pun and third 
Pjii order on the electrical field E: 
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X (^(ri,r) • ^*(r2,r)) E{T3,r)dTidT2dT3, (2) 

where x^^^ and e are the hnear electric susceptibility and the dielectric constant, 
x'^-* is the nonlinear susceptibility of third order, and we denote r = {x,y,z). 
We use the expression of the nonlinear polarization ([2]) , as we will investigate 
only linearly or only circularly polarized light and in addition we neglect the 
third harmonics term. The Maxwell's equations in this case becomes: 



IdB 

VxE^— — , 3 
c at 



V • D = 0, (5) 
V • 5 = V • i? = 0, (6) 

B = H, D = Pi,„+Pni, (7) 

where E and H are the electric and magnetic fields strengths, D and B are 
the electric and magnetic inductions. We should point out here that these 
equations are valid when the time duration of the optical pulses tg is greater 
than the characteristic response time of the media tq {Iq >> tq), and also when 
the time duration of the pulses is of the order of time response of the media 
{to < To)- Taking the curl of equation ^ and using Q and (O, we obtain: 



where A = is the Laplace operator. Equation ([5]) is derived without using 
the third Maxwell's equation. Using equation ^ and the expression for the 
linear and nonlinear polarizations ([T]) and ([2|) , we can estimate the second term 
in equation ^ for arbitrary localized vector function of the electrical field. It 
is not difficult to show that for localized functions in nonlinear media with and 
without dispersion W ■ E ^ and we can write equation ([8]) as follows: 
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We will now replace the electrical field in linear and nonlinear polarization on 
the right-hand side of ([9]) with it's Fourier integral: 



+00 



E{r,t)^ / E {r,Lu) exp {—iLut)dLj, 



(10) 



where with E (r, lo) we denote the time Fourier transform of the electrical field. 
We thus obtain: 



e{T — t)E (r, uj) exp {—iLOT)dujdT + 



\— 00 — 00 



(11) 
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X exp {—i {uj{ti — T2 + T3)))dLjdTidT2dT3. 

The causality principle imposes the following conditions on the response func- 
tions: 



e(r - i) = 0; x^^^ (n - t, ra - t, - t) = 0, 
T - t > 0; n - t > 0; i = 1,2,3. 



(12) 



That is why we can extend the upper integral boundary to infinity and use the 
standard Fourier transform 131: 



J £(r — t) exp {—iuJT)dT = J e{T — t) exp {—iuJT)dT 



(13) 
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The spectral representation of the linear optical susceptibility eo(tt') is connected 
to the non-stationary optical response function by the following Fourier trans- 
form: 
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e{(jj) exp {—iujt) — I e{T — t) exp {—i(jjr)dr. 



(15) 



The expression for the spectral representation of the non-stationary nonhnear 
optical susceptibility x^^"* is similar : 



hoo +00 +00 



x'-^Huj) exp {-icut) 



— 00 —00 —00 



X exp (-i {uj{ti - T2 + T3)))dTldT2dT3. 

Thus, after brief calculations, equation (fTT|) can be represented as 



(16) 



AE = 



E (r, Lu) exp {—iujt)dLO 



E{r,Lu) E {r,Lu) exp (—1 {ujt))duj. 



(17) 



We now define the square of the linear fc^ and the generalized nonlinear A;^; 
wave vectors, as well as the nonlinear refractive index 712 with the expressions: 



H^l — 2 — ft «2, 



where 



(18) 
(19) 

(20) 



The connection between the usual dimensionless nonlinear wave vector fc^; and 



the generalized one is: kf^i = k^i E (r, uu) 
equation ()17|) can be expressed by: 



In terms of these quantities. 



A£; = - / k'^{uj)E {r,uj)e-yiY>{~iLjjt)duj 
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(21) 



6 



Let us introduce here the amplitude function A{r, t) for the electrical field 



E {x, y, z,t) ^ A {x, y, z, t) exp {i{kQZ - wot)), 



(22) 



where loq and fco are the carrier frequency and the carrier wave number of the 
wave packet. The writing of the amplitude function in this form means that we 
consider propagation only in +z -direction and neglect the opposite one. Let us 

write here also the Fourier transform of the amplitude function A{r,Lo — cjo): 



+ CX3 



A[r,t) — / A{r^oj — uo)eyi'p{—i{ijj — LLiQ)t)(kij, 



(23) 



and the following relation between the Fourier transform of the electrical field 
and the Fourier transform of the amplitude function: 



E (r, uj) eiq){—iujt) = 
exp (— i (koz — uJat)) A{r,uj — ujq) exp {{{lu — LUo)t), 



(24) 



Since we investigate optical pulses, we assume that the amplitude function and 
its Fourier expression are time - and frequency- localized. Substituting p3|) .(|23 p 
and into equation (PT|) we finally obtain the following nonlinear integro- 
differential amplitude equation: 



AA{r, t) + 2ih 



(25) 



k'^{uj) ( 1 -|- n2(w) A{r,uj ~ /jJq) ) A{r,u} — ojq) exp (— — ujQ)t)duj 



Equation (|25p was derived with only one restriction, namely, that the amplitude 
function and its Fourier expression are localized functions. That is why, if 
we know the analytical expression of k'^iuj) and n2(w), the Fourier integral 
on the right- hand side of ()25p is a finite integral away from resonances. In 
this way we can also investigate optical pulses with time duration of the 
order of the optical period Tq = 27r/wo- Generally, using the nonlinear integro- 
differential amplitude equation (j25p we can also investigate wave packets with 
time duration of the order of the optical period, as well as wave packets with a 
large number of harmonics under the pulse. The nonlinear integro-differential 
amplitude equation (|25p can be written as a nonlinear differential equation for 
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the Fourier transform of the amplitude function A, after we apply the time 
Fourier transformation (l23l) to the left-hand side of (l25l) : 



AA (r, u! — uq) + 2iko 



OA (r, u! — ll>q) 
dz 



+ { [l + n2{uj) A{r,uj-iUQ) ] k^{Lu) ~ k^{iUQ)] A{r,uj - ujq) ^ 0. (26) 



We should note here the well-known fact that the Fourier component of the 
amplitude function in equation (|26p depends on the spectral difference Auj = 
ui ~ ojQ, rather than on the frequency, as is the case for the electrical field. 



3 From amplitude equation to the slowly vary- 
ing envelope approximation (SVEA) 

Equation (|25p is obtained without imposing any restrictions on the square of 
the linear fc^(w) and generalized nonlinear fc^^ = k'^ {uj)n2{uj) wave vectors. To 
obtain SVEA, we will restrict our investigation to the cases when it is possible to 
approximate k^ and fc^j as a power series with respect to the frequency difference 

Lo — loq as: 



fc2 (lo) = 



^'eo (^) ,2/ , , ^ (fe' (^o)) , s 
2 = k (ujo) H ^-r [uj - t^o) 



192 
2 



e (coo)) 



(uj - uJoY 



(27) 



i;M^^^^ii.M^'-^^i^-^.>^..^ (28) 

To obtain SVEA in second approximation to the linear dispersion and in first 
approximation to the nonlinear dispersion, we must cut off these series to the 
second derivative term for the linear wave vector and to the first derivative 
term for the nonlinear wave vector. This is possible only if the series ([27|) and 
(j28p are strongly convergent. Then, the main value in the Fourier integrals 
in equation ((25)) yields the first and second derivative terms in ((27)) . and the 
zero and first derivative terms in ((^5)) . The first term in ^7} cancels the last 
term on the left-hand side of equation ((^5)) . The convergence of the series 
((?7)) and ((23) for spectrally limited pulses propagating in the transparent UV 
and optical regions of solids materials, liquids and gases, depends mainly on 
the number of harmonics under the pulses .21j. For wave packets with more 
than 10 harmonics under the envelope, the series ((27)) is strongly convergent, 
and the third derivative term (third order of dispersion) is smaller than the 
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second derivative term (second order of dispersion) by three to four orders of 
magnitude for all materials. In this case we can cut the series to the second 
derivative term in ([7f|) . as the next terms in the series contribute very little to 
the Fourier integral in equation When there are 2 — 6 harmonics under 

the pulse, the series (P7)) is weakly convergent for solids and continue to be 
strongly convergent for gases. Then for solids we must take into account the 
dispersion terms of higher orders as small parameters. In the case of wave 
packets with only one or two harmonics under the envelope, propagating in 
solids, the series ([?f)) is divergent. This is the reason why the SVEA does 
not govern the dynamics of wave packets with time duration of the order of the 
optical period in solids. Substituting the series (P7)) and (^5)1 in (^5]) and bearing 
in mind the expressions for the time-derivative of the amplitude function, the 
SVEA of second order with respect to the linear dispersion and first order with 
respect to the nonlinear dispersion is expressed in the following form: 



Av4- 



(fcofc" + k'' 
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dt 



(29) 



where ko — fc(wo) and k^^i — k^i (wo). We will now define other important 
constants connected with the wave packets carrier frequency: linear wave vector 
fco = fc(wo) = ^o\/£{^o)/c; linear refractive index n{u!o) = y^eiujo); nonlinear 
refractive index n2{cJo) — 37rx^^''(^o)/£('^o); group velocity: 



t'(^o) = ^ = (30) 



2 y e d, 

nonlinear addition to the group velocity (feg^j 



P y _ 2fcon2 ^8712 

and dispersion of the group velocity fc"(wo) = d-^k/duj^^^^^. All these quantities 
allow a direct physical interpretation and we will therefore rewrite equation (|29p 
in a form consistent with these constants: 



OA 
'dt 



OA 



"2 



k^v dn2 



A 



A 



dt 



2fco 
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A kQvn2 



kov^ J dt'^ 



A 



A. 



(32) 
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This equation can be considered to be SVEA of second approximation with 
respect to the hnear dispersion and of first approximation to the nonlinear 
dispersion (nonhnear addition to the group velocity). It includes the effects 
of translation in z direction with group velocity v, self-steepening, diffraction, 
dispersion of second order and self-action terms. The equations ([5^ and (^9)) . 
with and without the self-steepening term, are derived in many books and papers 
[ni [m [H [m [H]. From equations ([321) , jlH) , after neglecting some of the 
differential terms and using a special " moving in time" coordinate system, it is 
not hard to obtain the well known spatio-temporal model. Our intention in this 
paper is another: before canceling some of the differential terms in SVEA ([5^ . 
we must write p2p in dimensionless form. Then we can estimate and neglect 
the small terms, depending on the media parameters, the carrier frequency and 
wave vector, and also on the different initial shape of the pulses. This approach 
we will apply in Section 5. As a result we will obtain equations quite different 
from the spatio-temporal ones in the femtosecond region. 



4 Propagation of optical pulses in vacuum and 
dispersionless media 

The theory of light envelopes is not restricted only to the cases of non-stationary 
optical (and magnetic) response. Even in vacuum, where e = 1 and Pni = 0, we 
can write an amplitude equation by applying solutions of the kind (|33p to the 
wave equation ([5]). We denote here by V{x,y, z,t) the amphtude function for 
the electrical field E{r,t) in vacuum: 



E (x, y, z,t) = V {x, y, z, t) exp {i{koz - ujot)), (33) 

where loq and fco again are the carrier frequency and the carrier wave number of 
the wave packet. We thus obtain the following linear equation for the amplitude 
envelope of the electrical field: 




' ''^ (34) 



The vacuum linear amplitude equation (VLAE) ([M]) is obtained directly from 
the wave equation without any restrictions. This is in contrast to the case of 
dispersive medium, where we use the series of the square of the wave vector and 
we require the series (|27l) to be strongly convergent. That is why equation (|34p 
describes both amplitudes with many harmonics under the pulse, and ampli- 
tudes with only one or a few harmonics under the envelope. It is obvious that 
the envelope V in equation (|34p will propagate with the speed of light c in vac- 
uum. Equation ([34]) is valid also for transparent media with stationary optical 
response e — const. In this case, the propagating constant will he v ~ c/y/sJI. 
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5 SVEA and VLAE in a normalized form 



Starting from Maxwell's equations for media with non-stationary linear and 
nonlinear response, we obtained an amplitude equation and a SVEA using only 
two restrictions, which are physically acceptable for ultra-short pulses. Having 
adopted the first restriction, namely, investigation of localized in time and space 
amplitude functions only, we introduced the amplitude equation (j25p . Following 
the second restriction, i.e., limiting ourselves with the case of a large number 
of harmonics under the localized envelopes, we obtained the SVEA ((32)) . As 
it was pointed out in the previous section, the second restriction do not affect 
the VLAE The next step is writing SVEA ^ and VLAE ^ in di- 

mensionless variables and estimating the influence of the different differential 
terms. In this case, the coefficients in front of the differential operators in (|32p 
and (|34p will be numbers of different orders, depending on the medium n and 
712, the spectral region of propagation fep and loq, the field intensity |Ao| , and 
the initial shape of the pulses, namely, light filament r± << zq (LF), Light 
Bullets (LB) r_L « zo or Light disks (LD) r± >> zq. With rx we denote here 
the initial transverse dimension, "the spot" of the pulse, and with zq we de- 
note the initial longitudinal dimension, which is simply the spatial analog of 
the initial time duration to, determined by the relation zo — vIq or zq = cto in 
the vacuum case. The SVEA ([32]) and VLAE ((34)) are written in a Cartesian 
laboratory coordinate system. To investigate the dynamics of optical pulses at 
long distances, it is convenient to rewrite these equations in a Galilean coordi- 
nate system, where the new reference frame moves with the group velocity for 
equation ((5^ . t' = t; z' = z — vt: 



I 



dA 
dt' 



n2 



kov dn2 



A 



A 



dt' 



1 



A 



A 



dz' 



d^A d^A 
2v- 



2fcn 



n2kov 



-A,A- 



2 2 V kov^ J \ dt'^ dt'dz' I 2 

and with the velocity of light for equation ((34)) . t' — t; z' — z — ct: 



(35) 



A 



A, 



With Aj_ 



* ~ 2fc^ ^ ~ 2k^cW^ dt'dz' 



(36) 



+ we denote the transverse Laplacian. We define the 
following dimensionless variables connected with the initial amplitude and with 
the spatial and temporal dimensions of the pulses through the relations: 



A^AqA''- V^VoV' 



r^x ; y = rj_y ; z 



zqz ; 
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t' = tor ; z = zqz" ; t = tof . (37) 



After the substitution of these variables in ([5^ , ([55)) , ([Ml) and ([55| and making 
use of the expressions for the diffraction Zdif = fegr^ and dispersion Zdisp = 
tg/k" lengths, we obtain the following five dimensionless parameters in front of 
the differential terms in the equations ([32| . ([35)) ([34)) and ([36)) : 



7 i;2 ^± o ^di/ J 2 2 I /I 1 2 



Zdi 



7 = A:o^o"2 l^oT 



71 = |Ao| ^2 + 



kov dn2 
2 



(38) 



Omitting the seconds in the new dimensionless variables and constants, the 
equations ([5^ . ([55)) . ([M)) and ([5S)) can be represented as follows: 
Case a. SVEA ([32)) in a laboratory frame (" Laboratory" ) 



/ 



2ia(5^ 



dA dA 



A 



dt 



Ai A + 



a^i d^A 



9z2 



/3 



at2 



+ 7 



(39) 



A. 



Case b. SVEA ([55)) in a frame moving with the group velocity: 



IdA . d^A 

V V // 

(40) 

Case c. VLAE (IMl) in a laboratory frame: 
Case d. VLAE ([55)) in a Galilean frame: 
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It should be noted here that equal dimensionless constants in front of the dif- 
ferential terms in both the "Laboratory" and "Galilean" frames are obtained. 
This gives us the possibility to investigate and estimate simultaneously the dif- 
ferent terms in the normalized equations (I39p , (HO]) , (311) and . We will now 
discuss these constants in detail, as they play a significant role in determining 
the diiTerent pulse propagation regimes. 

- The first constant a = fegZo = Sttzq/Ao determines with precision 2it the 
"number of harmonics" on a FWHM level of the pulses. Since we use the slowly 
varying amplitude approximation, a is always a large number (a >> 1). 

- The second constant = r\/ determines the relation between the initial 
transverse and longitudinal size of the optical pulses. This parameter distin- 
guishes the case of light filaments (LF) 5"^ = r\l Zq « 1 from the case of LB 
52 = r'i/z^ ^ 1 and the case of light disks LD 6'^ = r^/zg >> 1. For fight 
filaments << 1 and we can neglect the differential terms with coefficient 
(5^. It is not diflicult to see that in this case the SVEA (gO]) and VLAE 
PP]) . (pij) can be transformed to the standard paraxial approximation of the 
linear and nonlinear optics. If we set the possible values of the optical pulses 
transverse dimensions at 3 — Amm > r± > lOOfim, we can directly obtain the 
above distinction in dependence on the time duration of the pulses. For light 
pulses with time duration ns > to > AO — 50ps we obtain 6"^ = r'^/z^ « 1 
and we are in the regime of LF and paraxial approximation. In the case of 
fight pulses with duration to ~ 3 — Aps up to 500 — 600/s it is possible to reach 
(5^ = r'^/zQ = 1 and we are in the regime of LB. For pulses in the time range 
300/s — 30/s we can prepare the initial shape of the pulses to satisfy the rela- 
tion 5^ = t\Iz^ >> 1 and thus reach the LD regime. It is important to note 
here that the wave packets in the visible and UV ranges with time duration 
to > 30/s contain more than 10-15 optical harmonics under the pulse, so that 
we are still in SVEA approximation. In the last two cases (LB and LD) the 
differential terms with 6^ cannot be ignored and the equations p9)) . (liO]) and 
(PT|) governing the propagation of pulses with initial form of LB and LD are 
quite different from the paraxial approximation. 

-The third parameter is f3 = kor'^/zdis, where Zdis — ta/k" determines 
the relation between the diffraction and dispersion lengths. The dispersion 
parameter fc" in the visible and UV transparency region of dielectrics has values 
from k" ~ 10"^"'^ s'^/cm for gases and metal vapors up to fc" ~ 10^^^ s'^/cm for 
solid materials. It is convenient to express this parameter using the product of 
the second constant 6^ and the parameter (3i = kov'^k" by the relation (3 = PiS^. 
For typical values of the dispersion fc" in the visible and UV region listed above, 
the dimensionless parameter /3i is very small << 1), while for optical pulses 
propagating in the UV region in solids and liquids it may reach /3i oc 1. The 
parameter /3i can be also negative and may reach /3i ~ — 1 near electronic 
resonances and also near the Langmuir frequency in electronic plasmas. As it 
is shown in [17j . only in the case /3i ~ — 1 we can obtain the 3D+1 nonlinear 
Schrodinger equation from SVEA and ((40l) . 

- The fourth and fifth constants 7 — kQr\ri2 |^o|^ and aS'^^i ~ a5^n2 |^o|^ 



13 



are correspondingly the nonlinear coefficient and the coefficient of nonlinear 
addition to the group velocity (coefficient before the front of the first order 
nonlinear dispersion term). It is easy to estimate that for a >> 1 and 6^ > 1, 
we always have 7 >> a(5^7i. For optical pulses with power near the critical 
threshold for self- focusing 7 = 1 and less (linear regime) 7 << 1, the nonlinear 
addition to the group velocity is very small (ai5^7i << 1) and from here to the 
end of this paper we will neglect the terms with the first addition to the nonlinear 
dispersion. The analysis of the dimensionless constants performed above leads 
us to the following conclusion: Dynamics of wave packets with power near to 
critical for self-focusing 7 oc 1 in the visible and UV region in a media with 
dispersion are governed by the following SVEA equations: 
Case a. SVEA in laboratory frame (" Laboratory" ) 



Equations (j43ll and (|44l) are quite different from the well known paraxial spatio- 
temporal evolution equations. Here are included also the second derivative 
along the z direction, a mixed term and additional second derivative in time 
term. This leads to dynamics of the ultrashort fs pulses different from spatio- 
temporal model. In this paper we will investigate the propagation in linear 
regime, when 7 << 1. 



6 Fundamental solutions of the linear SVEA and 



The behavior of long pulses is similar to that of optical beams, since their 
propagation is governed by a equation where the nonparaxial terms become 
small. That is why we can expect the diffraction enlargement of long pulses 
to be of the same order as are the optical beams. The situation regarding LB 
and LD is different. Their propagation is governed by equations in media with 
non-stationary optical response - SVEA ([33]) and (gll), and by VLAE ([M]) and 
l|36p in media with linear stationary optical response (or vacuum), where the 
nonparaxial terms are of same order and bigger than transverse Laplacian. In 
this section we will solve the equations (|43l). (|44)) in linear regime (7 << 1) and 
will compare the solutions with the solutions of the linear VLAE (jM]) and ([36]) . 
Neglecting the small nonlinear terms in (|15]l . ([H]) we obtain: 




+ 7 A A. (43) 



Case b. SVEA in frame moving with group velocity ("Galilean"): 




(44) 



VLAE 
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a. Linear SVEA in a laboratory coordinate frame: 



b. Linear SVEA in a Galilean coordinate frame: 



For comparison we will rewrite here the corresponding linear VLAE: 
c. Linear VLAE in a laboratory frame: 



d. Linear VLAE in a Galilean frame: 



As can be expected, the equations for ultra-short optical pulses in vacuum 
and dispersionless media (|47l) . (fiSj) become identical with the equations with 
dispersion (I45p . (|46|) when the dispersion parameter (3i and ^^/Si are small. In 
equations (j45ll . (|46ll there are three dimensionless parameters, a, 6^ and /3i 
, while in ([T71) . (HS)) there are only two a and 5^. These parameters can be 
changed considerably in fs region, and this leads, as we can see later, to quite 
different dynamics in the particular cases. 

In the general case we apply the Fourier method to solve the linear SVEA 
([l5|) , ([35]) which describe the propagation of ultrashort optical pulses in medium 
with dispersion and linear VLAE (|T7)) . ([^5]) which govern the propagation of light 
pulses in vacuum and dispersionless media. We mark the Fourier transform of 
the amplitude functions of SVEA in Galilean frame (|45|) with Ao^kx, ky,kz,t) = 
F{A{x,y, z,t)), and in Lab frame ([46]) with AL{kx, ky, kz,t) = F{A{x,y,z,t)). 
The Fourier transform of the amplitude functions of VLAE in Galilean frame 
become (|T7)) . Bcikx, ky, kz,t) = F{V{x, y, z, t)), while in Laboratory frame (^5]) 
we write BL{kx, ky, kz,t) — F{V{x,y, z,t)). Applying spatial Fourier transfor- 
mation to the components of the amplitude vector functions A and V, the 
following ordinary linear differential equations in kx,ky, kz space for SVEA: 

a. Laboratory: 
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dt 

- (fc.' + ky^ + S^ih^ - 2aK)) Al - iPi + 1) (49) 
b. Galilean: 



r) A 



dt 

{kj + ky^ - S^Pikl) Ag - {f3i + 1) (50) 



and the following equations for VLAE: 
c. Laboratory: 



- ^laS'^ = {kj + k/ + S'{k^' - 2ak,)) 6'^, (51) 

d. Galilean: 



- 2^S' (a-k^)^^- {kj + ky') Bg - <5^^, (52) 

are obtained. We look for solutions of the kind of Al = A^^k^, ky, k^) exp(ir2it) 
and Ag = AG{kx,ky,kz)exp{i{lGt) for the equations ([SO)) , and for solu- 

tions Bl — B^ikx^ky, kz) exp(i$it) and Bg — Bcikx^ky, kz) exp{i^Gt) for the 
equations ([?T|) and ([5^ correspondingly. Let us denote the square of the sum 
of the wave vectors as: P = fc^^ + fcy^ + ^^(^z^ — 2akz)- The solutions exist 
when Ql, ^g, and $g satisfy the following quadratic equations: 



ni - 2 fir 



(53) 



^2^ _ g («-(/3l + l)fc.) ^^^ _ kx'+ky'-S^f3,Pz ^ Q 



(54) 



^2 r,/„. L ^^ _ kx + ky ^ ^ 



- 2(a - A:,)$G - 



<52 



(55) 
(56) 
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The solutions of ([551) . (EH) f^'' media with dispersion are: 



a 



,2 



/3i + l 



± 



(A + 1)' 



+ 



(57) 



Q - (/3l + 



± 



{a 



+ 1)2 <52(/?i + 1) 



(58) 



while the solutions of (|55p . and (|56p for dispersionless media and vacuum 
become: 



Now the necessity of a parallel investigation of the propagation of optical pulses 
in media with dispersion, in dispersionless media and in vacuum becomes obvi- 
ous. Further, we will introduce here the concept of weak and strong dispersion 
media depending on the value of dimensionless dispersion parameter /3i . When 
Pi << 1 we have media with weak dispersion. It is not difficult to calculate 
that in the optical transparency region of gases, liquids and solid materials, this 
parameter is usually very small ( /3i << 1). For media with weak dispersion the 
solutions of the characteristic equations with dispersion (fSTj) . (ISS]) are identical 
to the solutions for media without dispersion and vacuum ((55)) , (pD|) correspond- 
ingly. In media with strong dispersion Pi can reach the values /3i ~ 1 — 3 in 
the UV transparency region of solids and liquids. In this case, the solutions 
for media with dispersion will be slightly modified with the factor Pi + 1 with 
respect to the solutions without dispersion. We consider here the regime of 
propagation far away from electronic resonances and the Langmuir frequency 
in electronic plasmas, where it is possible to obtain a strongly negative disper- 
sion parameter Pi oc —1. We point out here again that in the case of LB, when 
Pi = —1, the amplitude equations ((5^ can be transformed into the 3D-f 1 linear 
and nonlinear vector Schrodinger equations jT7]. Generally said, the dispersion 
parameter Pi varies slowly from the visible to the UV transparency region of 
the materials from very small values up to /3i ~ 1 — 3 and that is why it does not 
influence radically the solutions and the propagation of optical pulses in linear 
regime. The other parameters a and change significantly. For example a 
varies from 10^ to 10^, while S'^ varies from 10"^ _ IQ^^ foj. Lp to 10° for LB 
and 10^ — 10"* for LD. This is the reason to investigate more precisely in the 
next paragraph the solutions of the equations for media with weak dispersion 
as air where Pi << 1, ([58]) . as we expect that the solutions for media with 
strong dispersion (UV transparency region of solids and liquids) will be only 
slightly modified by the factor Pi + 1; Pi < 1. We obtained solutions of the 




(59) 




(60) 



17 



characteristic equations ((57|) . ([5^ and ([SO)) . The solutions of the corre- 

sponding hnear differential equations SVEA (gH), ([501) and VLAE ([511), (ES]) in 
the fc-space become: 

a. Solution of SVEA in the k-space and laboratory coordinate frame: 



Al = AL{K,ky,k^,t = 0) X 



a 



b. Solution of SVEA in the k-space and Galilean coordinate frame: 

Ag = AG{k^,ky,k^,t = 0) X 

(62) 



, . / Q - (A + l)fc, / (a-(/3i + l)fc,)' fcg -t- fcg - J2/3ifc,^ 
c. Solution of VLAE in k-space and laboratory coordinate frame: 



Bl = BLik,,ky, k,, t = 0) exp \i ± ■\/a2 + ^7^2 j ^ j , (63) 
d. Solution of VLAE in the k-space and Galilean coordinate frame: 



Eg = BG{k^,ky,k„t = 0) exp f (a - k,) ± J + k^S^ t . (64) 



It is obvious that the solutions ([63| and ([64| of equations ([5T|) and (|52)) should 
be equal with accuracy - a wave number in z direction. This follows from the 
Fourier transform of such evolution equations and leads to only one difference 
between the solutions in the real space - the motion of the pulse in the z direction 
in a Laboratory frame and its stationarity in a Galilean frame. As it was 
pointed out at the beginning, we investigate here only localized in space and time 
initial functions of the amplitude envelopes. Thus, the images of these functions 
after Fourier transform in the kx,ky, k^ space are also localized functions. The 
solutions of our amplitude equations in k space (ISTI) . (|62p . (IS5)l and are the 
product of the initial localized in kx,ky, k-space functions and the new spectral 
kernels which are periodic (different from the Fresnel's one). The product of 
a localized function and a periodic function is also a function localized in the 
kx,ky, kz space. Therefore, the solutions of our amplitude equations in k space 
([M)) . (pT|) and are also localized functions in this space and we can 
apply the inverse Fourier transform to obtain again the fundamental localized 
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solutions in the x,y,z,t space. More precisely, we use the convolution theorem 
to present our fundamental solutions in the real space as a convolution of inverse 
Fourier transform of the initial pulse with the inverse Fourier transforms of the 
new spectral kernels: 

a. Fundamental solution of SVEA in laboratory coordinate frame: 



b. Fundamental solution of SVEA (l46ll in Galilean coordinate frame 




A{x,y,z',t') = F-^ (^AG{k^,ky,k„t = 0)^ 



c^-i I I . I « - + l)k, ^ /(a - + l)k,y kl + kl- 5-^l3ik1 , 



c. Fundamental solution of VLAE pT]) in laboratory coordinate frame: 

V{x,y,z,t)=F-^ [BL{k,,ky,k,,t^O)^ ® 



(^exp yi[a±\la'^ + kyS^j tj j . (67) 
d. Fundamental solution of VLAE in Galilean coordinate frame: 



V{x,y,z',t') = F-^ {BG{k^,ky,k,,t = Q)^ 

(^exp [i {{a -k,)± ^^^2 + P/JS^ , (68) 

where with we denote the spatial three-dimensional inverse Fourier trans- 
form and with ® we denote the convolution symbol. The difference between the 
Fresnel's integrals, describing propagation of optical beams and long pulses in 
linear regime, and the new integrals ([55)) . ([S7)) and ([55]) . which are solutions 
of the linear evolution equations (|l5)) . (|47p and ([^5]) is quite obvious. In 

addition, in the new spectral kernels there are three dimensionless parameters: 
a, and /3i. Let us fix a to be always large, i.e. a >> 1. As pointed out above, 
the condition a >> 1 is not necessary for VLAE, so we can investigate pulses 
with longitudinal duration of the order of the carrier wavelength in vacuum and 
in dispersionless media. To provide one qualitative analyze for the influence of 
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the other two parameters, and on the evolution of the initial pulse, we 
will rewrite the expression for the spectral kernel (jST]) of the solutions ([65]) of 
equation pS)) in the following form: 



/3i 



As a >> 1 and /3i < 1, the diffraction widening will be determined by the 
second term under the square root in ([55)1 : 

^ (fc' + (70) 



<52(/3i + 1) 

which determines the transverse diffraction and dispersion widening of the pulses. 
We pointed out above that the dispersion parameter varies very slowly within 
the limits < /3i < 10"'^, while the relations between the transverse and lon- 
gitudinal part varies sig nificantly 10""* < fp- < lO'^. This is why we estimate 
mainly the influence of the different values of on the diffraction widening. 
We investigate the following basic cases: 

a/ Long pulses, when «\. It is easy to estimate from (|70p . that the 
transverse enlargement k^'^ + ky'^ will dominate significantly as: 



1 



<52(/3i + l) 



» 1. (71) 



In the case of long pulses we have also aS'^ ~ 1 and as it can be seen from 
equations (j49[) and (|50[) we have similar (but not equal) diffraction length to 
that of optical beam (z^f^™ = kor^). The difference is only in the factors 
aS^ ^ 1 . When pulse propagate in optical transparency region of the air and 
gases f3i << 1, normalized dispersion parameter is to small that the main 
factor which determinate the diffraction widening become ^^"jy-'^ = ct'^^^di//' = 
fcgr^/zQ. The validity of this new diffraction formula for pulses will be studied 
more precisely in the next paragraph not only for long pulses but also for LB 
and LD. 

b/ LB: (5^ ~ 1. In the case of optical pulses with approximately equal trans- 
verse and longitudinal size, we obtain the following coefficient of the transverse 
kx^ + ky^ diffraction terms: 

^ 1/2. (72) 



s^Pi + 1; 



Hence, the diffraction and dispersion transverse enlargement will be reduced by 
the factor d'^{(3i + 1) with respect to the diffraction of long pulses and Fresnel's 
diffraction. 
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In addition, we will point out here an important asymptotic behavior of 
LB: When o? is small (pulses with only few harmonics under the envelope) 
and /3i << 1 (media with weak dispersion, dispersionless media and vacuum), 
the spectral kernels of the new equations tend to the asymptotical value ~ 

exp {i{'^k'^ + ky + [kz — aY)^^ = exp (i(|fc|t)), which is actually the spectral 

kernel of the 3D wave equation. For this reason we can expect for optical 
pulses with only one or two harmonics under the envelope (subfemtosecond and 
attosecond pulses) diffraction similar to the typical diffraction of the 3D wave 
equation, whose dynamics is characterized by internal and external fronts and 
a significant widening of the pulse. 

c/ Light disks: This is the case when the longitudinal size zq is mush shorter 
than the transverse size r± and >> 1. As indicated above, the typical time 
region for such pulses is 30 — 40/s < to < 200 — 300/s. We determine the 
lower limits of this relation from the condition >> 1, i.e a large number 
of harmonics under the envelope. This condition still holds true for pulses in 
the visible and UV regions with time duration 30 — 40/s. The dimensionless 
parameter in front of the transverse diffraction and dispersion (fca;^ + fcj,^) will 
be of the order of: 



We thus see that the transverse enlargement is of the order of 5^(/3i + 1), or 
negligible as compared with LB, and smaller by a factor of about ((5^(/3i + 
1))^ than in the cases of long pulses and Fresnel's diffraction. To summarize 
the results of this section, we can expect that the transverse diffraction and 
dispersion enlargement of LB should be smaller by a factor S'^{/3i + 1) than 
those of LP, while the transverse diffraction and dispersion enlargement of LD 
should be smaller by a factor of ((5^(/3i + 1))^ than those of long pulses and 
paraxial approximation. Practically no transverse enlargements of LD would be 
observed over long distances, namely, more than tens and hundred of diffraction 
length. 



7 Dynamics of light beam, LF, LB and LD in 
air 

In the beginning of this section we will discuss more widely the Galilean invari- 
ancy and connections between normalized equations written in different coordi- 
nate systems: Laboratory and Galilean. The Galilean invariancy of the SVEA 
([l5l) and VLAE (|T7)l is not obvious. Indeed, after using the transformation t' = t 
and z' = z — vt where v — 1 the new equations in Galilean frame, ([46]) and (|48p 
admit mixed 2, t terms and look quite different. The Galilean invariancy can be 
seen only from the kind of the fundamental solutions of equations (|65p - (|68p in 
Laboratory and Galilean frames. The solutions are equal with precision wave 
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z=0 



z=l 

ziir=Zdifl=7.S? cm. 



Figure 1: Intensity profiles of a Gaussian beam with initial condition 
Ax{x,y,z = 0) = exp^— ^-j^^ governed by the 2D paraxial equation ([7i|) . 

The transverse size (the spot) grows by factor \/2 over normalized distance 
z = 1. This correspond to a real distance zq = z'^ljY = '''•85 cm for the selected 
in the paper laser source on Aq — 800 nm. 

number , which gives the stationarity in Galilean and translation in z direction 
in Laboratory frame. The numerical solutions with initial conditions - Gaus- 
sian pulses provided in [22j for both coordinate systems demonstrate again that 
localized waves admit equal spatial and phase deformation in both coordinates 
and there is only one difference - stationarity in Galilean and translation with 
normalized velocity v=l in Laboratory frame. Naturally, the best way is to 
solve numerically the equations (|45l) or (|47l) in Laboratory system and to see 
one spatial and phase transformation of the pulse as well as its translation. In- 
convenience in such one approach is, that at long distances the pulses will move 
out of the grid. And here the Galilean invariancy of the normalized equations 
help us. It is well known that the times in Galilean and Laboratory systems are 
equal t' — t. As the normalized velocity is = 1, the same translation z = 1 in 
Laboratory frame correspond to time evolution in both system t' = t = 1 . This 
is demonstrated in [22] but here it gives us one additional opportunity. We can 
solve the equation in Galilean frame for long time t' without pulses to move 
out of the grid and to connect the normalized time t' in Galilean with normal- 
ized time t and translation z in Laboratory frame t' — t = z; v — \. After that, 
using normalized constants, we can obtain the real distance of propagation of 
optical pulses. It is not hard to see that normalized distance z — \ correspond 
to a real distance zq = fegr^ when a5'^ = 1. And this is the natural way to com- 
pare the dynamics of optical pulses governed by equation in Galilean frame 
with the evolution of a laser beam in scalar paraxial approximation described 
by normalized equation: 




(74) 
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where, as it is well known, z = 1 in normalized coordinates corresponds to a 
real distance zq = kQr\ called diffraction length. This length determines the 
distance where the laser beam increase its width on level from the maximum 
with factor-\/2. We investigate here only laser sources with spectrally limited, 
not phase modulated initial Gaussian profile. The optical lens and devices add 
additional phase modulation on the initial pulse and influenced on the widening 
in linear regime. 

Evolution of real laser pulse with the following characteristics is considered: 
light source form Tirsapphire with width on level e"^; r±_ = 100 ^m. Usually 
such small spot of the pulse is made by focusing by lens. To obtain no modulated 
in phase initial pulse the additional phase from the lens must be reduced to zero 
by a system of lens. The solutions of linear SVEA (|46p in Galilean frame are 
carried out for optical wave on wave-length Aq = 800 nm propagating in air and 
the following constants: carrying wave number fco = UbUj/c = 7.854 x lO'^; cm~^, 
where rif, « 1.0 for air; GVD coefficient fc" = 3.0 x 10~^^ sec^/cm; normalized 
GVD coefficient /?i = kgv'^k" = 2.1 x 10~^. To find the difference in dynamics of 
LF, LB, and LD we select different time duration of pulses for LF {to = 260ps), 
LB {to = 330/s) and LD {to = 33/s). Using the above parameters of the 
laser sources and the material constants we obtain the following dimmensionless 
parameters in SVEA (|46|) for the particular cases: 

a) long pulse and ai'^ = 1 (to = 260ps): a = 6.16 x 10^; 5^ ^ 1/6.16 x 10"^; 
/3i = 2.1 X 10-4. 

b) light bullet (330 fs): a = 785.0; = 1.0; /3i ^ 2.1 x 10"^ 

c) light disk (33 fs): a = 78, 5; = 100.0; /3i = 2.1 x 10"''. 

Other important parameter for comparing the pulse dynamics and paraxial 
evolution of a laser beam is the diffraction length Zdifj = koT^ = 7.85 cm. In 
addition, we should point out that all coming numerical computations are per- 
formed with pulses satisfying the boundary conditions fim^, ^^^^^ /2 ^(2;, y, z, t^x = 
and also lim^ j, ,2i-»±A/2 ^(fc^j ky^ kz,t) = 0, where L and A are respectively the 
spatial and wave-number intervals for the calculations. 

7.1 Evolution of Gaussian beam in paraxial approxima- 
tion 

The initial conditions for linearly polarized normalized Gaussian beam reads: 



yl — Ajr-X'^ 

A.(x,y,z = 0)=exp(^-^^±^^ . (75) 

The evolution of the initial Gaussian beam (j75p governed by the paraxial equa- 
tion (I74p is described by the Fresnel's integral or can be found by numerical 
calculation of the inverse Fourier transform of the solution in the {kx, fcj,)-space. 
The intensity profile of a solution A(x,y,z) of the paraxial equation (|74p with 
initial condition (j75p on a normalized distance z = 1 is illustrated on Fig.l. 
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Getting in ming the above real parameters of a laser system on 800 nm, the 
normalized distance z = 1 corresponds to one diffraction length and real dis- 
tance of zq = Zdiff = kor^^ = 7.85 cm. 




t'=fl 



=1: ZQ=2diff =7.-3? cm 



Figure 2: Transverse intensity distribution of 260 ps pulse on carrying frequency 
800 nm (long Gaussian pulse). Numerical solutions of the linear SVEA in 
Galilean frame is performed by the following particularly selected initial condi- 
tions to satisfy a6'^ = 1: Aj:{x,y,z,t = 0) = exp a = 6.16 x 10^; 

52 = 1/6.16x10-^ Pi = 2.1x10-''. The surfaces \A{x,y,z' = 0, i' = 0;t' = 
are plotted. The transverse size (the spot) grows by factor \/2 over normalized 
time-distance t' = z = 1. This correspond to real distance zq = z^lJY — 7.85 
cm equal to the diffraction length of a laser beam (compare with Fig.l). 



7.2 Evolution of long optical pulses (light filaments) 

The pointed above choice for the parameter aS^ = 1 of a long pulse is used 
particularly to compare it's diffraction with the diffraction length of a laser 
beam. We mark also that in the general case, the real diffraction length of a 
long pulse (ns or ps) is similar to the diffraction of laser beam and the difference 
is in the factor aS'^ or: 



pulse c2 beam i 2 A i 



(76) 



The validity of this expression is illustrated on the next two figures where the 
dynamics of an initial long pulse is governed by the linear SVEA (|^^ in Galilean 

-> ^ / 2 I 2 I /2 \ 

frame with initial condition A — A^x; Ax{x, y, z', t' = 0) = exp I — ^ j 

and dimensionless constants: 

a) long pulse on 260 ps: aS'^ = I; a = 6.16 x 10^ S'^ = 1/6.16 x 10"^ 
Pi = 2.1 X 10-*. 
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t'=C f=i t =S: 

zo=SzdHi<beainf=47.1 em 

Figure 3: Intensity distribution profiles of 43 ps pulse on carrying frequency 
800 nm (long Gaussian pulse). Numerical solutions of the linear SVEA (pS]) in 
Galilean frame is performed with initial conditions satisfy aS"^ = 6; a = 1. x 10^; 
(5^ = 6. X 10~^; (3i = 2.1 x 10^^. The transverse size (the spot) grows by factor 
1/2 over normalized time-distance t' = z — 6. This correspond to a real distance 

pulse _ c2 beam — A^J \ ™ 



b) long pulse on 43 ps: aS"^ = 6; a = 1. x 10 ^-^ — 6. x 10 ^; /3i = 
2.1 x 10-"*. 

Case a) is illustrated on Fig. 2, where the spot (x,y size) of the pulse is 
plotted and the parameters are selected to satisfy the relation aS^ — 1 . That is 
why the pulse enlarge its spatial width by a factor ^/2 on the same normalized 
time-distance t' = z = 1 as in the case of laser beam. Case b) is illustrated on 
Fig. 3., where the important dimmensionless parameter is aS^ — 6. As can be 
expected the spot of the pulse grows by factor ^/2 on six time longer distance 
than a laser beam and the special case a). Summarizing the results obtained 
for the linear regime of long pulses we conclude: The long optical pulses admit 
diffraction length of the same order to the one of a laser beam and this length 
can be equal only in some partial cases, satisfy aS^ = 1. 

7.3 Propagation of light bullets in linear regime 

The evolution of LB in media with dispersion, is governed by the same SVEA 
as in the case of long pulses. The shape of the LB is symmetric in the 
X, y and z plane, so that the linearly polarized initial Gaussian profile can be 
written as: 



A^x; a = 785; 5^ ^ = 1, /3i = 2.1 x 10"'' 

A^ix, y,z,t^O)^ exp (- + + V (77) 
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i'=785/2 



t'=785 
ZD=785 Zdifrtbea m) --62 m 



Figure 4: Evolution of a Gaussian light bullet with 330 fs time dura- 
tion governed by linear SVEA in Galilean frame under initial condition 

A^{x,y,z,t = 0) = exp(-^^-±2^±^), a = 785.0; S = 1.0; A = 2.1 x 10-*. 

The surfaces \A{x,y,z' = 0;i' = 0; t' = 785/2; t' = 785)^ are plotted. The 
transverse size (the spot) grows by factor \/2 over normalized time-distance 
t' = z = 785. For the selected in the paper laser source this corresponds to a 
real distance z^^lT = aS'^ z^'^'ij- ~ 62 m. 



From the qualitative analysis presented in the previous section, when 5^ = 1, 
the widening of the LB is expected to be a = 785; 2^""^^* = az^^'^J- = fc§rj[. 
The surface (x,y plane) of the solution of VLAE ([67| with initial conditions 
of the kind of ([77| on normalized time-distance t' — z — 785, calculated by 
exploit of FFT technique, is illustrated on Fig. 4. One can see that the pulse 
enlarge its spot by factor along a considerable distance of z^""/* = az'^'^^y = 
785z^f^™ = 6162.25 cm = 61 to. Let us remark once again that this result 
is only correct if the number of harmonics under the pulse, multiplied by 27r 
(dimensionless parameter a), is large. 

7.4 Dynamics of light disks in linear regime. 
Low-diffractive regime 

As it was mentioned in the beginning, optical pulses with small longitudinal 
and large transverse size, while at the same time the large number of harmonics 
under the pulse remaines, can be obtained without significant experimental 
difficulties. This can easily be realized in the optical region for pulses with 
time duration from 200 — 300 fs up to 30 — 40 fs. We consider again the 
propagation of LD in the framework of the solutions (pS)) of the SVEA ([^5]) in 
Galilean coordinates under initial conditions of the form: 
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t'=0 






^ Ml 




Xifi ~616 ED M 



Figure 5: Transverse intensity distribution of 33 fs pulse on carrying frequency 
800 nm (light disk) governed by the same SVEA and initial condition 
A4x,y,z,t = 0) = exp ^ ^ q, = 78.5; 6 = 100; f3i = 2.1 x 10-*. The 

surfaces \A{x,y,z' = 0;t' = 0; t' = 7850/2; = 7850)^ are presented. The 
LD enlarges its transverse size by factor over the normalized time-distance 
t' = z = 7850. This correspond to 7850 diffraction lengths of a laser beam or 
for the selected laser source: 2^4/^ ~ ^^"^^dtfT — ^^l m. 



A = A^x] a = 78.5; S"^ ^^-^ ^ 100, 

A^{x, y, z', t' = Q)= cxp + + ^ . (78) 

Results of calculations of solution with initial conditions of kind ([75]) , using 
FFT and inverse FFT technique are presented on Fig. 5. The numerical solution 
confirm our expectation that the LD enlarges its shape by factor y/2 on 7850 time 
longer distance than the diffraction length of a laser beam; 2^-|j = aS^z'^^jY = 
/cqT^/zq = 7850z^1^Y = 616 m. The new formula for diffraction length of optical 
pulses ([7S)) gives remarkable opportunity to select the parameters of the laser 
pulse and to obtain pulses with negligible diffraction. From ([75)) it is seen that 
z^ifY depends on the spot diameter of the pulse by four degree {z^lf'f ''D- 
If we use pulse large enough in transverse dimension we can obtain practically 
diffractionless pulses. For example, using ([76)) and light disk with waist r± = 1 
cm, time duration to = 33 fs (zo = 10 /^to), and ko = 7.85 x 10* cm~^ (Ao = 800 
nm) we can obtain pulse diffraction length of order of zflj'j: ~ 6160 km. Such 
LD will be propagate in transparency region of gases or vacuum on several 
thousand kilometers without practical diffraction enlargement. 
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8 Conclusion 



In this paper dynamics of ultrashort laser pulses in media with dispersion, dis- 
persionlcss media and vacniTim arc investigated in the frame of non-paraxial gen- 
eralization of the amplitude equation. In partial ease of media with dispersion, 
we obtained an integro - differential nonlinear equation, governing propagation 
of optical pulses with time duration of order of the optical period. The slowly 
varying envelope approximation (many harmonics under the pulse) reduced this 
amplitude integro - differential equation to the well known slowly- varying ampli- 
tude vector nonlinear differential equation with different orders of dispersion of 
the linear and nonlinear susceptibility. In case of propagation of optical pulses 
in dispersionless media and vacuum, we obtained an nonparaxial amplitude 
equation which is valid in both cases, namely, pulses with many harmonics and 
pulses with only one-two harmonics under the envelope. We normalized these 
amplitude equations and obtained five dimensionless parameters determining 
different linear and nonlinear regimes. The nonparaxial envelope equations for 
media with dispersion, dispersionless media and vacuum are solved in linear 
regime and new fundamental solutions, including the GVD, are found. In gases 
and vacuum the solutions of these equations predict new diffraction length for 
optical pulses ^^"j^*' = kgr^/zQ. We demonstrate by these analytical and numer- 
ical solutions a significant decreasing of the diffraction enlargement of fs pulses 
(LB and LD) in respect to paraxial widening of a laser beam and a possibility 
to reach diffraction- free regime. 
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